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Significance  points  for  some  tests  of  uniformity 
on  the  sphere 

J.  Keilson,  D.  Petrondas,  U.  Sumita,  J.  Wellner 

\ 


Abstract 


v/ 

v  Beran  (1968)  and  Gind  (1975)  have  proposed  several  omnibus  tests 
for  uniformity  on  the  unit  sphere  in  three  dimensional  Euclidean  space 
While  several  authors  have  contributed  to  providing  approximate  per¬ 
centage  points  for  the  limiting  distributions,  no  tables  of  the  limit¬ 
ing  distributions,  percentage  points  thereof,  or  finite  sample  dis¬ 


tributions  or  percentage  points  have  been  available.  In  this  paper 


•  -we-  fill  this  gap  by: 

(1)  finding  the  exact  distributions  of  the  statistics  of 
Beran  and  Gind  for  n  =  2; 

(2)  presenting  some  percentage  points  for  selected  small  and 
moderate  sample  sizes  obtained  by  Monte-Carlo  methods; 

(3)  evaluating  numerically  the  cumulative  distribution  func¬ 
tions  and  significance  points  of  the  limiting  distribu¬ 
tions  via  the  Laguerre  transform  method # (Keilson  and 

f'-  - 

Nunn  (1979),  Keilson,  Nunn  and  Sumita  (1981),  and  "'~'~ 
Sumita  (1981)). 
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§1.  Introduction 


Beran  (1968)  and  Gind  (1975)  have  proposed  several  omnibus  tests 
for  uniformity  on  the  unit  sphere  S  =  {x  e  R  :  |x|  =  1}  in  three- 
dimensional  Euclidean  space.  These  tests  are  consistent  against  all  alter¬ 
natives  and  are  locally  most  powerful  for  specific  alternatives.  Beran  and 
Gind  have  shown  that  the  limiting  distributions  of  these  statistics,  under 
the  null  hypothesis  of  uniformity,  are  those  of  weighted  sums  of  independent 
Chi-square  variables.  While  Prentice  (1975)  has  applied  the  methods  of 
Zolotarev  (1961)  and  Hoeffding  (1964)  to  provide  approximate  percentage 
points  for  the  limiting  distributions,  no  tables  of  the  limiting  distribu¬ 
tions,  percentage  points  thereof,  or  finite  sample  distributions  or  per¬ 
centage  points  have  been  available. 

Our  purpose  here  is  to  fill  this  gap  by: 

(1)  finding  the  exact  distributions  of  the  statistics  of  Beran  and 
Gind  for  n  =  2; 

(2)  presenting  some  percentage  points  for  selected  small  and  moderate 
sample  si:es  obtained  by  Monte-Carlo  methods; 

(5)  using  the  Laguerre  transform  method  -  Keilson  and  Nunn  (1979), 
Keilson,  Nunn  and  Sumita  (1981),  and  Sumita  (1981)  to  compute  the  cumula¬ 
tive  distribution  functions  and  significance  points  of  the  limiting  distri¬ 
butions  . 

In  Section  2,  we  summarize  the  result  on  the  limiting  distributions  due 
to  Beran  and  Gind  and  give  the  exact  distributions  for  n  =  2.  Section  3  con¬ 
tains  the  Monte-Carlo  results  for  finite  sample  sizes  and  description  of 
the  methods  used.  We  discuss,  in  Section  4,  the  numerical  procedure  for 
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evaluating  the  limiting  distributions  via  the  Laguerre  transform  method. 
An  application  is  given  in  Section  5  where  we  test  uniformity  of  orienta¬ 
tion  of  dendritic  fields  in  the  retinas  of  cats  subject  to  controlled 
visual  environments.  The  numerical  results  are  summarized  in  Section  6 
in  tables  and  graphs. 
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§2.  The  statistics:  limiting  distributions  and  exact  distributions  for 


Let  Xj,...,X  be  independent  and  identically  distributed  unit  vectors  in 

with  distribution  v  (so  v(S)  =  1,  where  S  =  {x  e  R^:  |x|  =  1}).  Let 

6. .  =  X.X.  =  arc  cos(X. *X.)  =  the  angle  between  X.  and  X.  for  i, j  =  1,2, .. .n. 

13  -l-j  -l  -y  6  -l  -3 

For  testing  the  null  hypothesis  that  v  is  the  uniform  distribution  on 
S,  Reran  (1968)  and  Gin£  (1975)  have  suggested  the  statistics 


Y°dd  ,  n'1  l  J  U  -  (2/ir)6  > 

i=l  j=l  J 


(2.1) 


n  n 


yeven  =  n’l  £  £  (j-  (2/ir)sin  6  } 


(2.2) 


i=l  j=l 


Y  =  yodd  +  yeven 
n  ~  n  n 


(2.3) 


Gin6  proposes  Y^ven  and  finds  that  Y°dd  is  simply  Beran's  form  of  Aine's  sta¬ 
tistic:  Y°dd  =  4T  (Beran,  1968).  It  is  also  known  that  Y°dd  =  4A0 

n  n  '  J  n  2,n 

(Prentice,  1978)  and  Yeven  =  G.  (Prentice). 

n  / 1  n 

OD 

Let  {2 j ) j _  j  be  independent  Chi-square  random  variables  with  2j  +  l  degrees 
of  freedom  (so  Z2k  -  Set 


a2k.l  *  C2k-l)'2[(j)  /k!]2 


(2.4) 


a ^  =  (2k- 1) " 1 (2k+2) ’ 1 [ (y)^/k ! ] 


(2.5) 


where  (y)k  =  y*  (y  ♦  1) - (y  +  k-1)  and  define 


y°dd  _  y  &2  j 

ki1  2k-lZ2k -1 


(2.6) 


Yeven  Y  2  . 

Y  a2kZ2k  ’ 


(2.7) 
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Y 


..odd  ..even 
\  +  \ 


(2.8) 


Theorem  1  (Gin£) :  If  v  is  the  uniform  distribution  on  S,  then 


lim  P[Y°dd  <  y]  =  P[Yodd  <  y]  , 

n-Ko  “ 

lim  P[Yn  £  y]  =  P[Y  <  y]  , 

n-«o 

and  lim  P[Y  <  y]  =  P[Y  <  v]  for  all  ye  R 

n-*=° 

Prentice  (1978)  uses  the  methods  of  Zolotarev  (1961)  and  Hoeffding  (1964) 

to  provide  approximate  percentage  points  for  these  limiting  distributions. 

Ke  will  compute  the  distributions  with  precision  via  the  Laguerre  transform 

method  [Keilson  and  Nunn  (1979),  Keilson,  N'unn  and  Sumita  (1981),  and 

Sumita  (1981)]  in  Section  4  and  present  tables  and  graphs  in  Section  6. 

Ke  now  consider  the  distributions  of  Y°dd,  yeven  and  Y  for  n  =  2. 

n  n  n 

Theorem  2:  If  v  is  the  uniform  distribution  on  S,  then 


P[Y°dd  <  y]  =  4(1  +  cos[y(y-2)]}  ,  0  <  y  <  2  , 

p^.even  <  yj  _  cos{arc  sin[y(l-y)]}  ,  1  -  7  <  y  £  1  , 

P[Y2  ^  )']  =  yd  +  cos  [g~ 1  Gy(3-y) )  ] }  ,  0  <  y  <  3  , 


-1 


where  g(t)  =  t  +  sint,  0  <  t  <  v,  and  g  denotes  the  inverse  of  g. 

Proof:  First  note  that  under  uniformity  P(01O  s  y]  =  (4n)~  /  /  sin6d?d4  = 

1  j  j  .  0  0 

y(l  -  cosy),  0  s  y  <  it.  Then,  writing  Y°  =  yU  ♦  1  ♦  (1  -  (2/ir)6j2^  + 

,.  2  .  -  r2sa  ..even  1,1  1  ,1  2  .  „  ,  ,x  2  .  .  . 

d  -  7  e21)>  =  2  -  (7)612,  Y2  =  2{2  +  2  +  (2  -  7  sln  6 12 3  +  (2  ‘  7  sin  021) 

2  2 

=  1  -  (— )sin  0^2*  a°d  Y2  =  3  -  (— )  ( © ^ 2  +  s*n  012^’  t*ie  statec*  distributions 
are  easily  obtained  by  straightforward  computation. 0 


For  future  reference,  we  record  some  moments  in  the  following  Proposi- 
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Proposition  3.  Under  uniformity  we  have 

E  0dd)  =  E(;Yodd)  =  i 

n 

E(  even  .  E(yeven  ,  1 
n  2 

E(Yn)  =  E (Y)  =  | 

while 

f.,odd-.  ,,  1.  16.  n-+o°  16.  ..  ,,  rv/odd, 

Var  (Y  )  =  (1  -  -)(2 - 7)  — >  (2 - j)  =  .37886  ...  =  Var[Y  ]  , 

7T~  7T 

Var(Y^Ven)  =  (1  -  ^-)  (— j  -  j)  ^  (— y  =  .04038...  =  VarfYeven] 

3tt“  3ir“ 

and 

VartY  )  =  (1  -  I)(I  -  -  ~)  =  .41924...  =  Var[Y]  . 

3tt  “  3tt" 

Furthermore,  the  third  and  fourth  cumulants  of  Y°dd  and  Yeven  are  given  by 
K,(Yodd)  =  .375219...  ,  K4(Y°dd)  =  .56252...  , 

K3(YCVen)  =  .0098016...  ,  K4(YeVen)  =  .00366374...  . 

Proof : 

The  finite  sample  means  and  variances  are  easily  obtained  by  elemen¬ 
tary  methods  upon  noting  that  0.^  and  0.^,  are  independent  for  j  *  j ' ;  the 
asymptotic  means  and  variances  follow  immediately  by  letting  n  -*■  ®.  The 
cumulants  of  Yodd  and  yeven  are  easily  computed  using  the  following  easily 
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§3.  Monte-Carlo  simulations 

In  our  simulation?,  a  function  subprogram  called  RAND  (University 
of  Rochester  file  #311*7*500,  Computer  Center)  and  the  IBM/360  computer 
are  used  to  generate  uniformly  distributed  random  numbers.  The  method 
of  Marsaglia  (1972)  is  then  employed  to  generate  points  from  a  distribu¬ 
tion  on  the  unit  3-sphere.  Using  n  such  random  points,  the  statistics 

yOdd  yeven  ancj  y  defined  in  Section  2  are  computed.  With  Monte-Carlo 
n  n  n 

samples  of  size  5000  for  n  =  100,  and  20000  for  n  =  5,  10,  20  and 

40,  the  percent  points  of  the  above  statistics  are  estimated  for  the 
significance  levels  a  =  .20,  .10,  .05,  .025,  .01  and  .001.  Those 
em^'^ical  finite  sample  percentage  points  are  further  smoothed  in  the 
following  manner.  For  each  level  a  and  sample  size  n,  let  a  smooth  func¬ 
tion  Y  be  defined  bv 
a  ,n 


Y 

a  ,n 


a  )  +  b  (-  - 
a a,n  n  2 
n 


(3.1) 


Here  Y  is  the  estimated  Y  value  for  each  n  at  level  a.  a  _  and  a 

a ,n  a, 2  a,® 

are  the  Y  values  at  level  a  for  n  =  2  and  n  =  °°,  respectively,  a  _  1S 

a,  i. 

found  from  Theorem  2,  and  a  from  the  Laguerre  transform  method  to  be 

’a,® 

described  in  Section  4.  b  is  the  estimated  slope  obtained  from  the 

a,n 

original  Monte-Carlo  results  via  the  straightforward  linear  regression. 

It  should  be  noted  that  the  smooth  function  (3.1)  coincides  with  the  known 
values  of  Y  when  n  =  2  or  n  =  ®. 

The  estimated  smoothed  percentage  points  of  the  three  statistics  are 
presented  in  Table  6.1  of  Section  6. 


§4.  Evaluation  of  the  limiting  distributions  via  the  Laguerre  transform 
Ke  have  seen  in  Section  2  that  Beran's  and  Gind's  limiting  statistics 
under  the  null  hypothesis  of  uniformity  are  infinite  sums  of  independent 

and  scaled  Chi-square  variates,  i.e.,  Yoc^,  Yeven  and  Y  as  given  in 
(2.6),  (2.7)  and  (2.8).  In  this  section,  we  discuss  numerical  evaluation 

of  Yodd,  Yeven  and  Y  and  propose  a  Laguerre  transform  approach  to  be 
described. 


r  2 

It  is  natural  to  decompose  the  infinite  sum  Y  =  [  aTZ.  in  (2.8) 

j  =  l  3  3 

into  two  parts,  the  sum  of  the  first  N  variates  and  the  remainder,  i.e. 


Y  =  S  +  V  •  S 
N  N  ’ 


N  00 

..  =  l  a2Z.  ;  VM  =  T  a2Z. 

N  ;  =  1  33  N  j-g+1  3  3 


(4.1) 


Correspondingly,  let 


,odd  _  ,-odd  yodd 
T  *  N  N 


yeven  _  seven  +  Veven 
N  N 


(4.2) 


where  S^t  and  Y^,  are  the  sums  of  odd  index  terms  of  and  V^,  respec- 

even  even  • 

tivelv.  The  variates  S,.  and  V.,  are  defined  similarly  for  even  index 

N  N 

terms.  S^.  is  the  finite  sum  of  independent  scaled  Chi-square  variates  and 
such  linear  combinations  may  be  regarded  as  positive-definite  quadratic 
forms  in  normal  variables.  Many  papers  have  been  published  on  the  distri¬ 
bution  of  such  quadratic  forms,  and  the  reader  is  referred  to  Johnson  and 
Kotz  (1970,  Ch.  29)  for  a  comprehensive  survey  of  the  literature.  An 
excellent  approach  to  the  numerical  evaluation  of  such  distributions  is 
that  of  Johnson,  Kotz  and  Boyd  (1967).  They  expand  the  distribution  func¬ 
tion  in  a  series  of  generalized  Laguerre  functions  and  evaluate  it  effi¬ 
ciently  by  taking  advantage  of  the  recurrence  relation  of  the  Laguerre 
functions.  When  this  procedure  is  applied  directly  to  the  distribution  of 
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S^,,  however,  one  encounters  numerical  difficulty.  The  coefficients 

■>  2 

a“  decreases  ranidlv  and  the  distributions  of  a  2.  become  very  concen- 

-1  J  J 

trated,  resulting  in  quite  slow  convergence  of  the  corresponding 

Laguerre  series.  Our  procedure  restructures  the  method  of  Johnson, 

Kot:  and  Boyd  and  overcomes  this  numerical  difficulty. 

The  Laguerre  transform  method  for  convolving  functions  has  been 

introduced  by  Keilson  and  Nunn  (1979),  Keilson,  N’unn  and  Sumita  (1981), 

and  further  studied  by  Sumita  (1981).  The  method  has  advantages  of 

accuracy  and  speed  which  make  it  an  attractive  candidate  for  problems 

of  this  type.  The  Laguerre  transform  method  has  peculiarities  and 

limitations,  however,  which  require  careful  refinement  for  particular 

contexts,  such  as  that  here.  The  reader  is  referred  to  three  basic 

papers  for  the  underlying  theory.  A  brief  summary  is  given  in  Appendix 

A  for  the  convenience  of  the  reader.  Our  basic  strategy  is  to  approxi- 

„  .  ..odd  ,  ..even  . 
mate  i  and  1  by 


vodd*  codd  ..odd*  Veven*  ceven  ..even*  , , 

Y\  =  SN  +  VN  ’  YN  =  S\'  +  VN  (4’3) 


odd*  ,  even* 

where  \^.  and  Y^,  are  the  Gamma  variates  having  the  same  first  two 
moments  of  and  V®ven,  respectively.  Correspondingly,  Y  is  approxi¬ 
mated  bv 


*  _  ,.oddv 
T\’  ~  N 


*  *rn*  *  \ 


(4.4) 


In  the  subsections  to  follow,  the  Laguerre  transform  procedure  is  des- 

odd*  PVPTl*  * 

cribed  for  evaluating  the  distribution  of  Y^  ,  Y^  and  Y^,.  The 
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validitv  of  the  results  is  also  examined. 


(A)  The  Laguerre  sharp  coefficients  of  the  Gamma  variate  r(q,  26) 
Let  F(a,  26),  a, 6  >  0,  be  the  Gamma  variate  with  p.d.f. 


g(x)  = 


a-1  26 

x  e 


0  <  x  < 


T (a) (26) 


(4.5) 


It  is  clear  that  the  variates  at 2^.  belong  to  this  family  with  a  =  j  +  y 

°  #00 
and  6  =  at.  Hence  the  Laguerre  sharp  coefficients  (g^)^  of  g(x)  pro¬ 
vides  a  basic  tool  for  the  procedure.  From  the  Laplace  transform 

oo  co 

Y(s)  =  /  e  SXg(x)dx  =  (1  +  26s)  a  and  the  identity  T#  (u)  =  £  g#un  = 

0  g  n=0  n 

Y(j  *  y7~) ,  those  coefficients  are  found  by  (cf.  Sumita  (1981),  Section 


g*  =  (l+6)~a  I  b 
ftn  n-n 

m=0 


(4.6) 


where 


b  =  n  (1 

n  i 

r=l 


~ )  .  n  >  1,  b0  =  1 


x  1  “  B  N  _  /•  .  l“Ctv  -  - 

Cn  '  W  ^  (1  -  >  n  >  1,  c0  =  1 

r=  1 


(4.7) 


The  accuracy  and  efficiency  of  the  Laguerre  transform  method  depend 
heavily  on  one’s  ability  to  represent  the  functions  present  with  a  sequence 
of  Laguerre  coefficients  of  reasonable  length  (say,  at  most  around  500 
coefficients  to  attain  5  digits  accuracy).  As  studied  theoretically  in 
Keilson,  Nunn  and  Sumita  (1981)  and  Sumita  (1981),  the  Laguerre  transform 
method,  when  applied  in  a  straightforward  manner,  cannot  tolerate  functions 
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too  closely  concentrated  at  zero  or  functions  too  great  in  extent. 

For  the  Gamma  variate,  this  point  can  be  observed  explicitly  in  (4.6) 

and  (4.7).  When  a  is  extremely  large,  b^  and  c^  become  so  large  in 

absolute  value  that  the  computer  may  not  tolerate  them.  The  other 

numerical  difficulty  arises  when  6  is  extremely  small  or  large.  In 

1-8 

this  case,  the  ratio  -r — -  becomes  verv  close  to  1,  and  one  would 

1  +  3 

expect  (cn)  to  have  a  long  tail.  This,  in  turn,  implies  a  long  tail 
of  g  .  Fortunately,  these  numerical  difficulties  can  be  avoided  by 
taking  advantage  of  the  divisibility  of  the  Gamma  variates  and  employ¬ 
ing  scaling.  In  brief,  the  first  difficulty  can  be  solved  through  the 
identity  y(s)  =  [ym(s)]M,  where  M  >  0  and  Y^(s)  =  (1  +  2Bs)  a//M  corre¬ 
sponding  to  tne  Gamma  variate  F(a/M,  2B) .  By  an  appropriate  choice  of 
positive  integer  M,  the  Laguerre  sharp  coefficients  of  T(a/M,  2B)  are 
obtained  with  reasonable  length.  Ke  then  convolve  them  M  times  on  the 
lattice  to  recover  the  Laguerre  sharp  coefficients  of  the  original 

Gamma  variate  T(a,  2S).  For  the  second  numerical  difficulty,  we  replace 

1  -c  8 

B  bv  cfi  so  that  the  ratio  becomes  well  below  1.  Then  (c  )  in 

1+cB  n 

(4.7)  decreases  rapidly  in  absolute  value  and  therefore  g^  decreases 
rapidly.  After  the  inversion  of  the  Laguerre  sharp  representation,  the 


proper  scale  factor  for  the  probability  density  function  is  restored. 

odd*  even*  * 

(B)  Algorithms  for  finding  the  sharp  coefficients  of  Y  ,  Y _ and  Y 


(1)  For  a  desired  small  variance  of  V^,,  select  N  for  S^,  and  V^.. 

Then  choose  the  dividing  factor  M  for  V^,  and  the  scale  factor  c  in  keeping 
with  the  conditions  of  (A).  The  number  L  of  Laguerre  coefficients  before 


truncation  will  be  discussed  in  (C)  below. 
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1  7 

(2)  Using  (4.6)  and  (4.7)  with  a  =  j  +  j  and  6  =  ca^ ,  calculate 
(g(.i ) *) q  oa“Z.  for  1  <  j  <  N. 

(3)  Obtain  (g°djJ#)o  of  cS°dd  by  convolving  (g(i)*)b  for  j  odd, 

.  ,  even#,L  ,,  „even  ...  , 

1  <  3  <  X.  Obtain  (g^,  n  )Q  of  cS^,  similarly. 

odd  ,  even  ,  ,  .  odd2  , 

(4)  Calculate  the  means  y^,  and  y^  and  the  variances  o^  and 

o-en2  of  VN°dd  and  V*ven ,  respectively,  from  (2.6),  (2.7)  and  Proposition 
3.  Using  (4.6)  and  (4.7)  with  a  =  (y°dd/o^d)  2/M  and  8  =  c*oVN  /2p\rN  ’ 
calculate  (h°dd* (M,N)  ) J.  By  convolving  (h°dd# (M,N)fl) J  M  times  with 


itself,  find  (h^#)o  °f  cVrU"-  0btain  ^^n'^O  °f  cVN^‘‘  similaTly- 

odd#  L  odd*  .  .  ,  odd#,L  ,  .,odd#,L 

(5)  Calculate  (f°  ||  )^  of  cY°a  by  convolving  (gN>n  )Q  and  (hN>n  )Q. 


even#, L 


.even*  ...  , 

similarly. 


Obtain  (f*,V*n#)J  of  cY*Ven*  similarly. 


#  L  *  .  ,_odd#,L  , 

(6)  Finally,  calculate  (f^  n)Q  of  cY^  by  convolving  (f^n  )Q  and 


..even#. L 
l\',n  V 

Remark 

ii  oo 

When  one  has  the  Laguerre  sharp  coefficients  Cf  )q  of  a  P*d.f.  f(x) 

ti  oo  t  , 

on  (0,ro),  the  inversion  of  (f  )q  to  the  values  of  f(x)  and  its  survival 

OO 

function  F(x)  =  /  f(x)dx  can  be  done  in  the  following  manner  (cf.  Keilson 
x 

and  Nunn  (1980)). 


(Rl)  Calculate  (^n(x))0  for  0  <  x  <  ®  by 

Wx)  *  sn-  rc2n-i-xiencxi  -  ^n-iwl  •  " s  1 

where  £^(x)  =  e  ^x. 

+  n  #  “  t 

(R2)  Obtain  f  =  [  f  and  calculate  f(x)  =  J  f  i  (x) 

n  Ln  m  n  n 


(R5)  Calculate  f*  =  -2  J  (-1)V+1+,  and  F(x)  =  j  ?^n(x)  . 

m=0  n=0 
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Ke  note  that  the  Laguerre  transform  bypasses  numerical  integration.  We 
also  note  that  the  algorithm  with  (Rl)  ,  (R2)  and  (R3)  produces  F^(x/c) 
where  F^(x)  =  F[Y^.  >  x]  . 

(C)  Validation  of  the  results 

For  the  calculation  of  the  distributions  needed,  the  scaling  factor 

is  taken  to  be  c  =  40.  The  two  values  N  =  2  and  N  =  4  are  used  and  the 

*  * 

dividing  factor  M  =  12  is  chosen  for  and  M  =  25  for  V..  The  length 

of  the  Laguerre  sharp  coefficients  is  L  =  502,  which  provides  12  digits 

2 

accuracy  of  the  p.d.f.  of  a.Z.  for  1  <  j  <  4. 

i  j  j 

There  are  two  different  factors  which  introduce  numerical  errors, 
the  truncation  of  the  Laguerre  sharp  coefficients  and  the  Gamma  approxi¬ 
mation  of  V  .  In  general,  it  is  quite  hard  to  quantify  truncation  error 
of  the  Laguerre  coefficients.  (Such  error  bounding  has  its  counterpart 
in  Fourier  series  theory,  where  error  bounding  is  known  to  be  extremely 
difficult.]  Theoretical  error  bounds  are  available,  so  far,  only  for  a 
certain  family  of  functions  (cf.  Keilson  and  Sumita  (1981)  and  Sumita 
(1981)).  Extensive  numerical  evidence,  however,  suggests  that  when  one 
chooses  L  large  enough  to  attain  a  given  accuracy  for  the  following  iden¬ 
tities,  then  the  function  values  are  likely  to  satisfy  the  same  accuracy. 


f(0+) 


-  u 
1  nf 

L  n 


(4.8) 


/  f(x)dx  =  ]  (-l)V 

0  n=0 


(4.9) 


/  xf(x)dx  =  4x  [  (-l)nnf 
0  n=0  n 


(4.10) 
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/  x2d(x)dx  =  16  x  l  (-l)nnV  .  (4.11) 

0  n=0  n 


#  SO9  * 

For  N  =  4,  (f^  n)p  ~  of  Y^  provides  10  digits  accuracy  for  all  equations 
(4 . S)  through  (4.10)  and  the  truncation  error  seems  to  be  negligible. 

Even  though  we  may  expect  the  Gamma  approximations  for  the  remainders 

V°dd,  Y-en,  and  V  to  introduce  little  error  since  Var(V^]  drops  rapidly 

-4 

(e.g.,  Var[V  ]  =  6.6  *  10  ) ,  no  analytical  justification  is  available, 

and  we  are  forced  to  take  indirect  means  for  testing  the  validity  of  the 

*  d 

approximation.  To  test  this  validity  we  note  that  Y„  -*■  Y  as  N  -*■  °°.  The 

•k  ★ 

c.d.f. 's  of  Y?  and  Y^  are  calculated  and  compared.  The  absolute  difference 

of  the  two  c.d.f. 's  is  found  numerically  to  be  bounded  by  1  *  10  5.  For 

a  second  accuracy  check  of  the  Gamma  approximation,  the  third  and  the 

odd*  even* 

fourth  cumulants  of  Y,  and  Y,  are  calculated  using  the  Laguerre 

4  4 

sharp  coefficients  and  compared  with  exact  values.  The  higher  moments 

of  Y°dd  and  Y^en  can  be  found  from  (4.10),  (4.11)  and  (R3)  of  (B) 

4  4 

#  f  4.  #  -f- 

with  the  identity  f  =  f  -  f  , ,  n  >  1,  and  f  =  f~.  The  cumulants  needed 

n  n  n- 1  n  0 

are  then  obtained  by 


K,  =  v. 


3u1u2  * 


K  =  p,  -  3p 


2u: 


2  -  4p1W3  + 


(4.12) 


12p1U2 


6p. 


where  p^  =  EfX1].  The  corresponding  exact  values  are  given  in  Proposition 
3.  The  computations  are  shown  in  Table  4.1.  The  absolute  difference  is 


bounded  by  3  *  10 
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Finally,  we  study  the  asymptotic  behavior  of  the  survival  function 
of  V,  and  compare  it  with  the  values  computed  via  the  Laguerre  transform 
method.  From  the  Laplace  transform  of  Y,  one  finds  through  asymptotic 
analysis  that 

-2x  X.  , 

Fv(x)  =  P[Y  >  x]  ~  - -  C 2Xn/2l  +  C—  +  A./f)  —  )  (4.13) 

Y  A  0  fi  1  G 

as  x  ->  +»,  where  XQ  w  4.95221  and  X}  "  -4.27140.  Details  are  given  in 

Appendix  B.  In  Figure  4.2,  the  asymptotic  expansion  of  Fy(x)  in  (4.13) 

* 

is  plotted  with  the  survival  function  of  Y^  derived  via  the  Laguerre 
transform  method.  Their  absolute  difference  is  found  to  be  bounded  by 
1  x  10  ^  for  x  >  3.5. 

A  similar  approximation  for  the  remainder  of  an  infinite  sum  of 
independent  variates  is  employed  in  Sumita  (1979)  to  evaluate  the 
multiple  convolutions  of  the  Logistic  variates.  There  a  direct  accuracy 
check  is  possible  and  accuracy  to  seven  decimal  places  is  attained. 

All  calculations  were  carried  out  on  a  DEC10  computer,  in  a  time¬ 
sharing  mode  using  APL  as  the  programming  language.  Relevant  formulae 
are  coded  in  a  straightforward  manner,  with  no  attempt  made  to  optimize 
the  subroutines  for  speed  and  accuracy.  In  spite  of  this,  the  results 
displayed  here  were  obtained  with  CPU  times  in  seconds  with  no  evidence 
of  numerical  problems. 
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K3  (Y0dd) 

K3  (Y°dd*) 

0.3752193777 

0.3752164153 

K4  (Yodd) 

K4  fY°dd*) 

0.5625201543 

0.5625200535 

Absolute 

Difference 

0.0000029624 

Absolute 

Difference 

0.0000001008 

K3  (Yeven) 

K3  (Y®ven*) 

0.0098016768 

0.0098007718 

K4  (Yeven) 

K4  (Yfen‘) 

0.0036637420 

0.0036637195 

Absolute 

Difference 

0.0000009050 

Absolute 

Difference 

0.0000000225 

Table  4.1. 

The  cumulants:  Yodd 

and  Y  vs. 

vodd*  .  Veven* 

_ aJ.d  Y.i 

Figure  4.2 
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§5 .  Application 

In  an  experiment  to  determine  the  effect  of  different  visual  stimuli, 

six  cats  were  divided  into  three  groups  of  two,  and  each  group  was  subjected 

to  a  different  visual  stimulus:  horizontally  polarized  light  (H) ,  vertically 

polarized  light  (V) ,  and  unpolarized  or  "normal"  light  (N) .  Orientations  of 

the  dendritic  fields  were  then  measured  at  15  to  16  sites  in  the  retinas  of 

each  of  the  six  cats.  The  complete  data  set,  in  coordinates  described  by 

Figure  5.1  below,  is  given  in  Appendix  C. 

The  orientation  of  the  dendritic  fields  of  the  two  cats  exposed  to 

normal  light  is,  presumably,  uniform,  with  no  preferred  orientation,  and 

the  question  is:  what  is  the  effect  of  polarized  light?  Hence  we  wish  to 

test  the  null  hypothesis  of  uniformity  of  orientation  (on  the  unit  hemisphere) 

of  the  dendritic  fields  of  the  H  and  V  groups.  To  test  this  hypothesis,  we 

«  n  n  .  *  a  a 

use  Gina's  statistic  yeven  =  n"  £  £(■=■-  (— )  sin  X.X.),  where  X.X.  is 

n  i=l  j=l  1  1  3 

the  angle  between  the  observations  X^  and  X^ .  Since  X^  is  represented  by  a 

vector  x^  =  [cos$^  sin0^,  sin<f^  sin0^,  cos0^]  with  Ijx^l^  =  1,  one  easily 

A  y 

finds  that  cos  X.X.  =  x.x.  and  therefore 

l  j  -i-j 


(5.1)  X^Xj  =  arc  cos[jcos(0^  -  6  ^ )  { 1  +  cos(<f>^  -  <f>j)} 
+  i-cos(0i  +  0j){l  -  cosOfK  -  6j)}]  . 


7 

/ 


I 


"elevation" 


Figure  5.1 


A  summary 

table 

of  the 

P-values  using 

Table  6.1 

is 

given  below. 

Cats 

n 

Yeven 

n 

P- 

value 

Cats 

n 

Yeven 

n 

P-value 

H1 

15 

.4508 

>> 

.20 

H1  +  H2 

31 

.9181 

M  .040 

H2 

16 

.7321 

.120 

V1 

16 

.4923 

>> 

.20 

V1  *  V2 

32 

.8603 

*  .056 

V2 

16 

.6239 

> 

.20 

N1 

15 

.6260 

> 

.20 

N!  +N2 

31 

.4421 

»  .20 

16 

.6076 

> 

.20 

Table  S.l 

As  the  table  shows,  the  test  of  uniformity  is  not  significant  for  any  indi 
vidual  cat.  Grouping  Hj  and  Hj  or  and  V^,  however,  does  give  significant 
results,  whereas  grouping  Nj  and  Nj  is  still  not  significant. 


In  Table  6.1,  the  estimated  smoothed  percentage  points  for  various 
values  of  n  and  significance  levels  a  are  presented,  using  the  method 
described  in  Section  3.  For  n  =  2  and  n  =  °°,  the  values  are  taken  from 
the  formulae  in  Theorem  2  and  Table  6.2,  respectively.  For  comparison 
purposes,  we  record  the  approximate  significance  points  of  Prentice 
(1978).  The  quantiles  of  the  limiting  distributions  of  Y0C^,  Yeven 
and  Y  are  given  in  Table  6.2,  where  a  =  P[Y  >  y] ,  etc.  More  detailed 
values  of  the  distributions  are  presented  in  Table  6.3,  with  step  size 
0.05.  Figure  6.1  and  Figure  6.2  show  the  graphs  of  the  c.d.f.  and  the 
p.d.f.  of  Yoc^,  yeven  and  Y.  All  the  values  associated  with  the  limit¬ 
ing  distributions  are  calculated  by  the  Laguerre  transform  method  des¬ 
cribed  in  Section  4. 
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Table  6.1 


Selected 

Smoothed  Percentage 

_  .  .  .  ..odd  ..even 

Points  for  Y  ,  Y 

n  n 

and  Y 

n 

Sample 

Significance 

level  a 

Si:e 

Statistics 

.20 

.10 

.05 

.025 

.01 

.001 

Y°dd 

n 

1.4097 

1.5903 

1.7129 

1.7978 

1.8725 

1.9597 

2 

yeven 

0.6180 

0.7225 

0.8012 

0.8585 

0.9102 

0.9715 

(Exact) 

n 

Y 

n 

1.9004 

2.2084 

2.4359 

2.5991 

2.7458 

2.9195 

yOdd 

n 

1.4171 

1.7788 

2.0771 

2.5456 

2.6704 

3.3946 

S 

yeven 

n 

0.6269 

0.7345 

0.8348 

0.9323 

1.0598 

1.3838 

Y 

n 

1.9132 

2.2767 

2.6154 

2.9387 

3.3338 

4.3277 

yOdd 

n 

1.4162 

1.8063 

2.1534 

2.4802 

2.8914 

3.8606 

10 

yeven 

n 

0.6353 

0.7490 

0.8567 

0.9619 

1.1000 

1.4476 

Y 

n 

1.9278 

2.3124 

2.6803 

3.0389 

3.4927 

4.6366 

yOdd 

n 

1.4151 

1.8135 

2.1832 

2.5385 

2.9935 

4.0913 

20 

yeven 

n 

0.6405 

0.7582 

0.8696 

0.9777 

1.1183 

1.4657 

Y 

n 

1.9372 

2.3327 

2.7137 

3.0866 

3.5652 

4.7610 

yOdd 

n 

1.4144 

1.8155 

2.1959 

2.5654 

3.0424 

4.2060 

40 

yeven 

n 

0.6434 

0.7634 

0.8765 

0.9858 

1.1270 

1.4712 

Y 

n 

1.9423 

2.3434 

2.7306 

3.1098 

3.5997 

4.8156 

yOdd 

n 

1.4140 

1.8161 

2.2029 

2.5808 

3.0712 

4.2747 

100 

yeven 

n 

0.6452 

0.7666 

0.8809 

0.9908 

1.1320 

1.4735 

Y 

n 

1.9455 

2.3501 

2.7409 

3.1236 

3.6199 

4.8460 

yOdd 

1.41563 

1.81631 

2.20727 

2.59079 

3.09000 

4.32040 

00 

yeven 

0.64643 

0.76879 

0.88384 

0.99413 

1.13534 

1.47452 

Y 

1.94776 

2.35459 

2.74772 

3.13268 

3.63309 

4.86522 

Prentice 

(1978) 

(♦4) 

yOdd 

- 

1.7750 

2.1750 

- 

3.0500 

4.3000 

(*16) 

Yeveji 

- 

0.7625 

0.8833 

1.1313 

1.4680 
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Table  6.2.  The  Quantiles  of  the  Limiting  Distributions 


a 

Y  (  O  t*  I*  ) 

V (EVEN) 

Y 

0.999 

0.18563 

0.16005 

0.47887 

0.995 

0.22343 

0.18567 

0.55003 

0.990 

0.24414 

0.20184 

0.58984 

0.975 

0.28309 

0.22580 

0.65944 

0.950 

0.32650 

0.24992 

0.73038 

0.900 

0.39123 

0.28234 

0.82748 

0.850 

0.44743 

0.30846 

0.90397 

0.800 

0.50073 

0.33171 

0.97212 

0.750 

0.55346 

0.35351 

1.03633 

0.700 

0.60686 

0.37460 

1.09897 

0.650 

0.66181 

0.39553 

1.16151 

0.600 

0.71904 

0.41674 

1.22509 

0.550 

0.77939 

0.43850 

1.29091 

0.500 

0,84376 

0.46113 

1.35992 

0.450 

0.91324 

0.48501 

1.43340 

0.400 

0.98921 

0.51060 

1.51289 

0.350 

1,07352 

0.53849 

1.60028 

0.300 

1,16916 

0.56940 

1.69861 

0.250 

1.28015 

0.60472 

1.81215 

0.200 

1.41363 

0.64643 

1.94776 

0.150 

1.58259 

0.69823 

2.11887 

0.100 

1.81631 

0.76879 

2.35459 

0.050 

2,20727 

0.88384 

2.74772 

0,025 

2.59079 

0.99413 

3.13268 

0.010 

3.09000 

1.13534 

3.63309 

0.005 

3.46337 

1.23944 

4.00709 

0.001 

4.32040 

1.47452 

4.86522 

(e.g.,  P [Y  >  0.47887]  »  0.999) 


ci  ci  ci  ci  ri  ci  c*  Ti  ci  ci 
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Table  6.3.  The  Cumulative  Distribution  Functions  of  Yot^,  yeven  and  y 


o.oo 
o.os 
0.10 
0.15 
0.20 
0.25 
0.30 
0.35 
0.40 
0.45 
0.50 
0.55 
0.60 
0.65 
0.70 
0.75 
0.80 
0.85 
O.VO 
0.95 
1.00 
1.05 
1.10 
1.15 
1.20 
1.25 
1 .30 
1.35 
1.40 
1.45 

1.50 
1.55 
1.60 
1.65 
l./O 
1.75 
1  .HO 
1  *U5 
1  .VO 
1.V5 

.00 
.05 
.  10 
.15 
.20 
•  25 
.30 

.40 

.45 

2.50 


W -OPb ( X ) 

W -EVC H ( K ) 

r<K) 

K 

F-0»6(K) 

r.KvcH(K) 

r  <:<> 

0.00000 

0.00000 

0.00000 

2.55 

0.97308 

1.00000 

0.92B99 

0.00000 

0.00000 

0.00000 

2.60 

0.97542 

1.00000 

0.93500 

0.00000 

0.00000 

0.00000 

2.65 

0.97756 

1.00000 

0.94051 

0.00005 

0.00025 

0.00000 

2.70 

0.97952 

1.00000 

0.94557 

0.00180 

0.00887 

0.00000 

2.75 

0.98131 

1.00000 

0  •  95021 

0.01147 

O.OSOOB 

0.00000 

2. BO 

0.9G293 

i.ooooo 

0.95446 

0.03353 

0.13223 

0.00000 

2.85 

0.98444 

1.00000 

0.95836 

0.06670 

0.24167 

0.00000 

2.90 

0.98581 

1.00000 

0.96194 

0.10739 

0.36068 

0.00006 

2.95 

0.98706 

1.00000 

0.96522 

0.15235 

0.47607 

0.00040 

3.00 
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APPENDIX  A 

The  Laguerre  Transform 


The  Laguerre  polynomials  Ln(x),  defined  by  the  Rodrigues  formula 
L  (x)  *  ex(^)n(xne‘x) ,  form  a  set  of  orthonormal  polynomials  with 
weighting  function  w(x)  =  e  x  on  (0,  »)  (see,  e.g.,  Szego  (1975)).  The 
associated  Laguerre  functions  £n(x)  =  e  L^Cx)  then  provide  an  orthonormal 
basis  in  L2(0,  •) .  For  any  f(x)  e  L2(0,  ») ,  one  has  the  Fourier-Laguerre 
expansion 

CD  00 

(A.l)  f(x)  =  l  fit  (x)  if  l  =  f  f(xK  (x)dx  . 
n=0  n  n  n  0 

to  oo  dcf  00 

Let  T^(u)  =  l  f+un  and  T*(u)  *  l  f*un  I  (I-u)tJ(u).  Since  £  ^(x) u11 
*  n=0  *  n=0  n=0 

(l-u)-1exp(-^x(l+u) (1-u)"1},  one  has,  when  f(x)  is  integrable  on  (0,  ®) , 


00  00  oo 

t!(u)  -  (1-U)  j  fV  =  (1-u)  l  un  /  f(x)£  (x)dx 
*  n=0  n  n=0  0  n 


*  /  f  (x)exp{-Jjx(l+u)  (l-u)_1}dx  ,  0  £  u  <  1  , 

0 

i  .e. , 


(A. 2)  tJ(u)  • 


where  <t(s)  «•  /  e'sxf(x)dx  is  the  Laplace  transform  of  f(x).  Let  f(x)*g(x)  * 

0 

x 

/  f(x-_y)g(y)dy.  Since  (s)  *  $f(x)$  (s) , 

0  1  8  1  8 
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(A.3)  T*,g(u)  =  T*(u)T*(u) 

for  functions  f(x),  g(x)  that  are  both  in  1^(0,  «)  and  L2(0,  ®) .  From  (A. 3) 
one  obtains 


(A. 4)  (f*g)#  =  l  f#  g# 

s  n  L .  n-msm 
m=0 


The  transformation  via  (A. 3)  maps  functions  f  (x)  ,  g(x)  into  sequences 

(fn)(gn)»  and  their  continuum  convolution  f(x)*g(x)  is  mapped  into  a  lattice 

t  n  # 

convolution  and  then  back  onto  the  continuum  via  fT  =  7  f  and  the  repre- 

n  Ln  m  1 

m=0 

sentation  (A.l).  This  transformation  procedure  was  introduced  originally 
in  Keilson  and  Nunn  (1980).  The  Laguerre  transform  was  extended  subse¬ 
quently  in  Keilson,  Nunn  and  Sumita  (1981)  to  handle  functions  on  the  full 


continuum,  and  further  studied  by  Sumita  (1981). 
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APPENDIX  B 

Asymptotic  Behavior  of  Fy(x)  =  P[Y  >  x]  as  x  -*•  +°° 

In  this  appendix,  we  derive  the  asymptotic  expression  of  Fy(x)  as 
x  -*■  +”.  The  Laplace  transform  of  the  p.d.f.  of  Y  is  given  by 


(B.l) 


*  (s,  -  n  :  9 

j=l  j  J  2a^ 


Let 


(B.2) 


“<s)  -  (epH)372  i 

1  J=2  J 


Ke  note  that  A(s)  is  regular  for  Re(s)  >  -07.  Let  Fy(x)  =  P[Y  >  x]  and. 


e.x  «  -(w  -  e.)x 

define  C'(w)  =  L{e  F..(x)}  -  /  e 

0 

one  has 


Fy(x)dx.  Then,  from  (B.l)  and  (B.2) 


1  -  <?Y(w  -  6  )  1  -  a(w  -  6.) 

(B.5)  V(w)  =  - - —  =  - - —  •  A(w  -  0 1 )  ♦ 


w  -  e 


l 


l 


v 


l  -  x(w  -  e}) 

w  -  0, 


After  a  little  algebra,  one  finds,  from  (B.2),  that 


(B.4) 


i  -  a (w  -  e^)  /e^ 


w  -  6 . 


.3/2 


n-1  n-3 

“  n  1.  2  2 

I  (-D  (?) 


On  the  other  hand,  A (w  -  6^)  and 


/6jW  n=3 

1  -  A(w  -  ex) 
w  -  8~ 


are  regular  for  Re(w)  >  ■  e2 


with  6j  -  02  <  0.  Hence,  one  has  the  Taylor  expansion 


A(w  -  0.)  *  [  A  w71  . 

n=0 


(B.S) 
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